Sub-threshold production of E(1385) baryons 
in Al+Al colhsions at 1.9A GeV 
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First measurement of sub-thresliold E(1385) production is presented. Experimental data are 
presented for Al+Al reactions at 1.9A GeV measured with the FOPI detector at SIS/GSI. The 
E(1385)/A ratio is found to be in good agreement with the transport and statistical model pre- 
dictions. The results allow for a better understanding of sub-threshold strangeness production and 
strangeness exchange reaction which is the dominant process for K~ production below and close-to 
threshold. 

PACS numbers: 25.75.-q, 25.75.Dw 



Relativistic heavy ion collisions at SIS energies provide 
interesting opportunities for studying hot and dense nu- 
clear matter. It allows to address fundamental aspects of 
nuclear physics such as the nuclear-matter equation-of- 
state [l|,ll, H, 01 and the question whether hadron prop- 
erties undergo modifications in such an environment [5|. 

Indications for in-medium modifications of charged 
kaon production and propagation have been experimen- 
tally observed by the KaoS and the FOPI 0, i] col- 
laborations at SIS energies. This beam energy range (1 
- 2A GeV) is particularly well suited for studying in- 
medium properties of strange particles since, as they are 
produced below or close-to threshold, their production 
process is sensitive to nuclear in-medium effects. While 
sub-threshold are mostly produced via multi-step 
processes, transport models indicate that sub-threshold 



K production results from strangeness exchange reac- 
tions TT + Y^iir^+B whose rate is intimately linked to 
the hyperon (Y = A, S) yield 

In addition, recent calculations based on the chiral 
theory predict an important coupling of the K~ to 
the medium via S(1385), A(1405) and A(1520) reso- 
nances 

[m, [H H El and reveal in particular, the 



importance of the E(1385) hyperon in sub-threshold 
production TB^I. Theoretical predictions on the in- 
medium cross-section of this strangeness exchange reac- 
tion differ widely [HI, [Hill Q and, on the other hand, 
experimental data on E(1385) are scarce. 
We report in this letter on the first measurement of 
charged S(1385) (called S* in the following) production 
in Al+Al collisions at a beam kinetic energy of GeV. 
This measurement corresponds to sub-threshold produc- 
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tion since the threshold beam kinetic energy for S* pro- 
duction in elementary reaction is -Ethr. = 2.33 GeV. The 
measured yield is compared to the predictions of statis- 
tical and transport models. The only other available 
measurement of E* production yield in heavy ion col- 
lisions has been reported recently by the STAR collab- 
oration for the Au-|-Au system at the top RHIC energy 
{^/JWn = 200 GeV) [i3,[l3l- 

The experiment has been performed with the FOPI 
detector at the SIS Heavy-Ion Synchrotron of GSI- 
Darmstadt by using an Al beam of kinetic energy of 
1.9 A GeV on an Al target. The beam intensity was 
chosen to be 8 • 10^ ions/s and the target thickness was 
567 mg/cm^. The FOPI detector is an azimuthally sym- 
metric apparatus made of several sub-detectors which 
provide charge and/or mass identification over nearly the 
full solid angle. The central part of the detector set- 
up, placed in a super-conducting solenoid, consists of a 
Central Drift Chamber (CDC) surrounded by a barrel 
of plastic scintillators. The forward part is composed 
of two walls of plastic scintillators and a second drift 
chamber. For the present analysis, A and S* were re- 
constructed from their decay products measured in the 
CDC (23° < diab < 114°). These are identified by 
their mass, which is determined from the correlation be- 
tween magnetic rigidity and specific energy loss. More 
details about the configuration and performances of the 
different components of the FOPI apparatus can be found 
in 1^, 1^, 2^ . The events were centrality selected by im- 
posing conditions on the multiplicity of charged particles 
detected in the forward part of the detector. The results 
presented in the following correspond to the most central 
collisions {ageo < 315 nib) representing about 20% of the 
total geometrical cross section. 

In order to assess the P(i;*- -I- E*+)/P(A + 1]°) ratio, 
yields of E*^, as well as primary A have to be de- 
termined. Note that the reconstructed A yield includes 
decay A from E*^ which cannot be isolated with the FOPI 
detector because E° decays into a A and a photon. While 
primary A are reconstructed from the invariant mass 
of (p,TT^) pair candidates (see [2l|, for more details 
about A reconstruction in FOPI), E* are reconstructed 
from a topological analysis of their double two body de- 
cay: 



A 



P 



(BR = 88%, cr = 5 fm) 
(BR = 64%,cr = 7.89 cm). 



Since the short life-time of the E* does not allow to dis- 
antengle its decay vertex from the collision vertex, its 
reconstruction consists in correlating A with primary pi- 
ous. Therefore, in order to minimize the combinatorial 
background in the (A, tt^) pair candidates, a high purity 
sample of A is extracted from the data. These A are re- 
constructed from identified pions and protons with the 
following conditions: 

i) a cut on momentum {piab > 0.1 GeV/c) is applied 



in order to reject particles which are spiraling inside the 
CDC; 

ii) the distance of closest approach (DCA) between 
tracks and the primary vertex in the transverse plane 
(DCAp > 0.8 cm, DCA^- > 1.9 cm) is used to enhance 
the fraction of protons and pions coming from a sec- 
ondary vertex; 

iii) a cut on transverse flight distance {dt < 20 cm) and 
on the pointing angle of the A (A0 < 4°), the latter being 
the difference between the geometrical and kinematical 
azimuthal angles; 

iv) a cut on transverse momentum of the A in the lab- 
oratory {pt > 0.3 GeV/c) is used to increase the purity 
of the signal. 

The invariant mass of the particle pair is calculated 
from its four-momenta at the intersection point of the 
two tracks. The corresponding distribution is shown in 
the inset of Fig. [TJ A total of about 10^ A are recon- 
structed with a signal-to-background ratio close to 10 
from the analysis of 290 million events. E*^ are then 
reconstructed by calculating the invariant mass of these 
A with charged pions after applying a two-sigma mass 
selection cut around the A nominal mass (shown by the 
dashed lines in the inset of Fig. [T]). 
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FIG. 1: Invariant mass spectra of pn^ pairs (inset, upper 
panel) and Avr^ pairs. The solid histogram and crosses de- 
note the data and the scaled mixed-event background, re- 
spectively (upper panel). The lower panel shows the signal 
after background subtraction. The following characteristics 
of the signal are shown: number of counts in the signal (S), 
signal-to-background ratio (S/B) and significance (SIGNIF). 
The parameters extracted from the fit to the data (mean mass 
value (MEAN) and the width (F)) are also reported. 



The reconstructed E*~ and E*+ are shown in Fig. [TJ 
The combinatorial background is obtained with the 
event- mixing method [23|. For that purpose, the two 
decay particles (A and tt^) are taken from two differ- 
ent events which present the same particle multiplicity 
in the CDC. In addition, the two events are aligned to 
the reaction plane in order to have the same reference 
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system for both particles. The reaction plane is esti- 
mated event by event utilizing the transverse momentum 
procedure detailed in [2J]. In the invariant mass range 
1.35 to 2 GeV /c^, 16 excited states of S's have already 
been observed [25|. Therefore, the normalization of the 
mixed background was done in the mass range 1.25 - 
1.32 GeV/c^ (grey zone in the upper panel of Fig. [1]) 
in order not to bias the reproduction of the combinato- 
rial background. The shape of the resulting mixed-event 
background describes the combinatorial background and 
is indicated by crosses in Fig. [1] (upper panel) . The verti- 
cal bars of the crosses correspond to the statistical errors. 
After background subtraction (Fig. [1] lower panel) , the 
remaining peak in the mass spectrum is fitted with a 
Breit-Wigner function. Within the width interval (F), 
about 3100 S*^ are found for the applied set of cuts 
in the analysis. The mean mass value and width ex- 
tracted from the fit are in a good agreement, within sta- 
tistical errors, with the values reported by the Particle 
Data Group [25|. The width of the E* is expected to 
increase with baryonic density [5^ but the accuracy of 
our measurement doesn't permit to draw any conclusion 
about a possible broadening. On the right side of the 
T,*^ signal, some excited states are visible as it has been 
observed in the invariant mass distribution presented by 
the STAR collaboration [il,[i3|- 

The losses due to decay, acceptance and reconstruc- 
tion efficiency have to be corrected in order to extract 
the particle yields. These corrections are determined by 
means of extensive GEANT simulations modeling the 
full detector response. The IQMD model 27| is used 
as generator for the underlying heavy- ion event. The 
E*^ and A resonances are embedded in those events 
(one signal per event) and they are generated separately 
with a momentum distribution according to the Siemens- 
Rasmussen formula (28j which describes an expanding 
system with a temperature T and a radial expansion ve- 
locity /?. The choice of the values of these parameters 
(T = 90 MeV, /3 = and 0.3) is imposed by previ- 
ous measurements 2^ 3^. Afterwards, a full simulation 
of the detector, including resolutions in energy deposi- 
tion and spatial position, Front-End-Electronics process- 
ing and tracking, is performed and the resulting output 
is subject to the same reconstruction procedure as for the 
experimental data. Simulated and experimental spectra 
of all relevant quantities were carefully compared. Since 
no significant differences were found, we conclude that 
the apparatus is properly described by the simulation 



311 . 132l |. Finally, the reconstruction efficiency is deter- 
mined by computing the ratio of reconstructed particles 
in the simulation to those initially embedded into the 
background events. The systematic error on the A and 
E*^ yields was evaluated in two steps. First the recon- 
struction efficiency was estimated from simulated A and 
E*^ with two values of the radial flow velocity {(3 = 
and /? = 0.3). Then the A and E*^ signals were recon- 



structed under different sets of conditions on the relevant 
quantities previously discussed (j>iab, DC A, dt, Ac/) and 
Pt). In order to minimize the statistical errors, the yields 
of E*+ and E*~ are summed and the following ratio can 
be extracted: 



P(E* 



E*+) 



P(A + EO) 



= 0.125 ± 0.026(stat.) ± 0.033(syst.). 



Note that using such a ratio has the additional advantage 
that the influence of the chosen momentum distribution 
for the simulated signal is cancelled out to some extent. It 
is worth to point out that the measured ratio vary, within 
errors, by less than a factor two from SIS to RHIC despite 
the different processes involved in the E* production [16|, 

a- 

Our experimental ratio is compared to a statistical 
model predictions for Al-I-Al collisions The results 

are shown in Fig. [2] where other particle species such as 
proton, pion, A and K*^ are included. 
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FIG. 2: Measured hadron yield ratios compared to the ther- 
mal model calculations for Al-I-Al collisions at 1.9A GeV. The 
symbols are data and the lines are the model calculations. 

The statistical model, based on the canonical ensem- 
ble, reproduces the ratios with a temperature of 76 MeV 
and a baryonic chemical potential of /is = 816 MeV. This 
chemical temperature is lower than the kinetical temper- 
ature of 95 ± 10 MeV derived from previous measure- 
ments [2^, H^l- The same model, used to extract ther- 
mal parameters from Au-|-Au collisions at ^/snn = 2.7 
GeV reproduced particle ratios with T = 64 MeV 
and /is = 760 MeV. It is surprising to notice that at the 
same energy, a lighter system (Al-I-Al) brings out higher 
temperature and baryonic chemical potential, while the 
volume represents 10% of the one of the Au-|-Au system. 
Finally, a higher temperature appears better suited to de- 
scribe particle ratios when including a baryonic resonance 
and further studies on other strange particles including 
mesonic resonances as (j) and K{892) will be useful to con- 
strain the temperature. As it has been reported in (34] . 
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an additional parameter as 75 which takes into account a 
non fully chemical equilibrium of strange particles could 
also be used to describe all the particle ratios. 

The measurement is also compared to the prediction 
from the UrQMD transport model 35| for Al+Al colli- 
sions [11]. It is worth to point out that this model is 
the only transport code predicting E* production at SIS 
energies. The results of comparisons with statistical and 
transport models are summarized in Tab. IH 



Yield ratio 


Data 


Therm. Model 


UrQMD 


P(E*"+E* + ) 
P(A+E") 


0.125 ± 0.042 


0.097 


0.177 



TABLE I; Experimental yield ratio and predictions from sta- 
tistical and transport models. The error of the experimental 
result corresponds to the quadratic sum of statistical and sys- 
tematic errors. 



be reached at SIS energies. 
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The transport model prediction for Al-I-Al collisions 
is in relatively good agreement with the measurement 
when taking into account the upper limit of the error 
bars. The dominant process involved in the UrQMD 
model in order to produce S* is the fusion of A and 
pions (76%) with an average cross section of about 37 
mb [36|. The other processes used in the transport 
code to produce E* are the interactions between E and 
pions (12%) and, N* and A resonances with baryons 
(12%). The time dependence of these reactions exhibits 
a maximum around 8 fm/c which is in agreement with 
the predictions of the IQMD transport model [l3| on the 



time production of strange particles such as kaons 37|. 
During the time span between the E* production and 
the chemical freeze-out, about 8% of pions coming from 
the decay of the E* are lost in inelastic interactions. 
This marginal loss could be attributed to the small size 
of the system created in Al-I-Al collisions. 



In summary, we have presented new results about the 
P(E*- + E*+)/P(A + E") ratio in Al+Al collisions at 
1.9^ GeV. For the first time, the E*^ resonances were 
measured 400 MeV below their production threshold. 
The experimental result was compared to the predictions 
of a statistical model and a transport model. The mea- 
surement is in good agreement with the transport model 
prediction where the dominant reaction to produce the 
E* is the fusion of A and pions. On the other hand, the 
comparison with the thermal model suggests that a chem- 
ical freeze-out temperature of about 80 MeV is needed to 
describe the measured ratio. 

The reported measurement should be used as an 
input for other transport codes for testing calculations 
on sub-threshold strangeness production. In addition, a 
measurement of the system size dependence of the E* 
production, as well as other strange resonances, will give 
the opportunity to assess effects of a partial restoration 
of the chiral symmetry at high baryon density that can 
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